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C/^ I Abstract 

■ In the present paper we are studying scale properties of twist-2 conformal operators in 

supersymmetric Wess-Zumino model. In particular, we are interested in a construction of mul- 
I tiplicatively renormalized conformal operators. We show, that in order to find multiplicatively 

renormalized operators in this model, it is sufficient to find multiplicatively renormalized oper- 
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(NJ ■ ators only for one member of operator supermultiplet. We found a closed analytical solution 

I for fermionic conformal operator, which together with supersymmetry transformations could be 

' used to find remained multiplicatively renormalized bosonic operators. Moreover, we found, that 

■ the knowledge of fermionic diagonal and non-diagonal anomalous dimensions matrices allows 

i us completely reconstruct the forward anomalous dimensions matrix in singlet case. 

iu'. 1 Introduction 

^ . Parton distributions in QCD satisfy the Balitsky-Fadin-Kuraev-Lipatov (BFKL) pp and Dokshitzer- 
• ^ ! Gribov-Lipatov-Altarelli-Parizi j2] equations. When considered in supersymmetric limit these equa- 
^ I tions reveal some remarkable properties. For example, in the maximally supersymmetric TV = 4 
d ' Yang-Mills theory there is a deep relation between BFKL and DGLAP evolution equations In 
particular, the anomalous dimensions of Wilson twist-2 operators in = 4 SYM could be found 
from the eigenvalues of the kernel of BFKL equation. Moreover, the anomalous dimensions matrix of 
the twist-2 operators in A^ = 4 Supersymmetric Yang-Mills theory has a very simple form with its el- 
ements being expressed in terms of one universal anomalous dimension with shifted arguments Oil]. 
In leading order, universal anomalous dimension was found to be proportional to the ^'-function, 
what may be considered as an argument in favor of integrability of this model jS]. Generally, the 
conformal invariance of the theory in leading order allows us to construct multiplicatively renormal- 
ized quasi-partonic operators [H] up to this order in perturbation theory. However, in next-to-leading 
order in perturbation theory multiplicative renormalization of conformal operators is violated due to 
necessity in regularization of arising ultraviolet divergences. 

The Efremov-Radyushkin-Brodsky-Lepage (ER-BL) equation [Zj may be viewed as some kind of 
the generalization of the DGLAP equation for the case of non-forward distribution functions. In this 
case in and out hadronic states in matrix elements of conformal operators are different, what allows 
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us to study scale properties of hadron wave functions. In the latter case the matrix elements of corre- 
sponding conformal operators are considered between vacuum and hadronic state. It is known, that 
the eigenf unctions of ER-BL evolution equation are directly related to multiplicatively renormalized 
conformal operators. Moreover, knowing ER-BL evolution kernels one may find both diagonal and 
non-diagonal parts of anomalous dimensions matrix for corresponding conformal operators. 

Recently, there was a great interest in properties of before mentioned operators in the context 
of famous AdS/CFT correspondence 0. Namely, there are some calculations of the anomalous 
dimensions of such operators in the limit of large j (Lorentz spin) from both sides of the AdS/CFT 
correspondence [31121 EE]- There are, also some predictions for anomalous dimensions of other 
types of operators from string theory J2] , partially confirmed by field theory calculations. It should 
be noted, that up to this moment only diagonal part of anomalous dimension matrix have been 
studied in the context of AdS/CFT correspondence and it would be interesting to compare the 
results for its non-diagonal part with the appropriate result from string theory. So, we believe, that 
non-diagonal parts of anomalous dimensions matrices in supersymmetric theories may reveal some 
interesting properties and thus are worth to study. 

In this paper we start our study with a simple supersymmetric model - supersymmetric Wess- 
Zumino model ^3] which contains only matter superfields. This model allow us to understand general 
properties of the ER-BL equations and their solutions in the supersymmetric theories. As the next 
step we suppose to perform a similar analysis for field theories with more supersymmetries. We 
will be mainly interested in the construction of multiplicatively renormalized conformal operators 
for non-forward kinematics and subsequent understanding of their string counterparts. This will be 
a subject of one of our next papers. 

To be precise, in the present work we calculate evolution kernels for non-singlet and singlet 
unpolarized conformal operators together with similar quantities in polarized case. We study the 
transformation properties of twist-2 conformal operators in Wess-Zumino model under supersymme- 
try transformations. At this point, it is natural to introduce fermionic operators (with respect to 
quantum numbers), as components of corresponding supersymmetric chiral multiplets. It is shown, 
that the latter allow us reduce the problem of finding the solution of evolution equations in singlet 
case (here we have operator mixing adding extra complexity) to the similar problem for fermionic 
operators. The corresponding problem for fermionic operators turns out to be much simpler and 
we will be able to obtain a closed analytical solution for multiplicatively renormalized fermionic 
operators, which together with supersymmetry transformation could used to obtain a solution of 
evolution equations in the singlet case. Moreover, we found, that the knowledge of fermionic diago- 
nal and non-diagonal anomalous dimensions matrices allows us completely reconstruct the forward 
anomalous dimensions matrix in singlet case. 

Our present consideration is based on two approaches developed for studying non-forward partonic 
distribution functions. The first one follows the articles of Mikhailov and Radyushkin [Til 1^ IT^- 
where an effective method for calculation of the evolution kernels of ER-BL equation was proposed 
US] . The properties of evolution kernels were later studied in [12] . The absence of the gauge 
fields in our problem considerably simplifies the whole calculation, so here, we are following Ref. [T^ . 
where the scalar 0^g^ theory in the space-time dimension six was considered. The second approach is 
based on papers of Miiller ^3021; where the renormalized conformal Ward identities were extensively 
exploited to find non-diagonal parts of anomalous dimension. 

The present paper is organized as follows. In section 2, after briefly recalling some necessary 
background about distribution amplitudes, conformal operators and supersymmetric Wess-Zumino 
model, we give our results for ER-BL evolution kernels in unpolarized and polarized cases together 
with some comments on the details of this calculation. In section 3 we introduce an approach for the 
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evaluation of non-diagonal parts of anomalous dimensions matrices for conformal operators, we are 
interested in, based on the analysis of conformal Ward identities for Green functions with operator 
insertions. The obtained non-diagonal anomalous dimensions matrices were checked to coincide with 
those coming from ER-BL evolution kernels evaluated before. In section 4, we combine our conformal 
operators into one supersymmetric chiral multiplet and study the constrains on anomalous dimensions 
of conformal operators following from supersymmetric Ward identity. Next, in section 5 we construct 
multiplicatively renormalized fermionic operator up to next-to-leading order in perturbation theory 
and explain how corresponding solutions could be found for other operators. Finally, section 6 
contains our conclusion. 



2 ER-BL evolution kernels in the supersymmetric Wess- 
Zumino model 

2.1 Preliminary 

The field-theoretical background for the study of the distribution amplitudes is provided by their 
relation to the matrix elements of non-local operators sandwiched between hadron and vacuum states 

m = 7^ / dz^e"P^^~ < 0|^(0)^(^_)|/i(p+) >, (1) 



27r „ 

where p+ is hadron momentum, x - momentum fraction carried by field (p and z_ stands for light-cone 
coordinate. 

ER-BL equation describes the scale dependence of the distribution amplitudes (DA) or more gen- 
eral non-forward parton distribution in contrast to scale dependence of forward distribution functions 
governed by DGLAP equation. The ER-BL equation is generally written as 

^^) = j dyY{x,y\a{ii))(j){y, /i) , (2) 

where evolution kernel 

V(x,y|a) = £v(°)(a;,y) + V«(a;,y) + ... (3) 

is given as a series in the coupling constant and is calculated in perturbation theory. In leading order 
the solution of evolution equation © could be constructed introducing moments of DA with respect 
to Gegenbauer polynomials in the case of bosonic (on quantum numbers) DA or Jacobi polynomials 
in the case of fermionic DA. For example, for non-singlet bosonic DA we have: 

'rfxCf (x-x)0(x,Q2) = ml,{^^')\PY;tQ., (4) 



where x = 1 - x and {Q\Ot^^{^l'^)\PY^'^ = {0\Ol^{^'^)\P) / P^^^ denotes reduced expectation value of 
non-singlet bosonic conformal twist-2 operator: 

Ot = l^^A^d^Yl^cf (^^) (5) 

where T> = d — d, d = d + d. In leading order this and other conformal operators, we will consider 
later, are multiplicatively renormalized, so that the solution of evolution equation can be easily 



3 



found. The evolution equation for DA Q leads to renormalization group equations for corresponding 
conformal operators: 



^ k=0 



where the anomalous dimensions matrix 



/ a \ (0)^ , / a \ 2 m 



7..= (^j7lX.+ (^j + (7) 

is diagonal in leading order and could be extracted from the expression for evolution kernel: 

/■I j 

/ dxCf^{x - x)V{x, y\a) = ^ 7jfc(a)C^^^(y - y). (8) 

So, we see that the problem of determining the evolution kernels for DA and the problem of 
calculation of anomalous dimensions matrix for corresponding conformal operators are closely re- 
lated. In what follows, we are using two different approaches to determine DA evolution kernels and 
anomalous dimensions matrices of our conformal operators. Equations, similar to (jH)) will be used 
only to check, that both methods give us the same results. 



2.2 Wess-Zumino model 

The field content of supersymmetric Wess-Zumino model ^H] consists of the (anti)chiral supermulti- 
plets containing off mass shell a complex scalar field ip, one Majorana fermion field ip ("quark" field) 
and a complex auxiliary scalar field F: 

^y,e) = ^{y) + V29ij{y) + 99J^{y) (9) 
<!>{x,9,9) = ^{x) + V29ij{x) + 99J^{x) 

- i9a^9d^,^{x) - -99998 ^d^'^ix) - ^999a^,d^'i){x) (10) 

The lagrangian for this model is conveniently written using superfield formalism and is given by 

C = j d^9dP9^,^' J d^9Yijk^%^<!?^ j d¥F^^'^l>i<lj<lfe , (11) 

where Yijk is Yukawa coupling, symmetric with respect to all indices. The first term in the la- 
grangian stands for kinetic terms of all fields present in this model, second and third terms denote 
superpotential (and its hermitian conjugate), which remains non-renormalized to all orders in pertur- 
bative theory. Due to non-renormalization of the superpotential the /3-function for Yukawa coupling 
completely determined by the wave function renormalization of fields. To simplify notation it is 
convenient to put Y^^'' = g. The first coefficient of Yukawa /5-function /3((7^) = boj^^ + . . . in 
Wess-Zumino model is bo = —3 and the leading order anomalous dimension of the chiral superfield 



is 74. = 7^2 



(47r)2 

In what follows, we will also need this lagrangian written in component notation. Rewriting 
complex scalar field ip = {A + iB) j \/2 in terms of real scalar A and pseudoscalar B fields, complex 
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auxiliary scalar field = (F + iG) / \/2 in terms of real scalar F and pseudoscalar G auxiliary fields 
the Wess-Zumino lagrangian takes the form 

+ g[F{A^ + B^) + 2GAB-ii){A + -i^B)ij]. (12) 

The supersymmetric transformations for the component fields could be summarized as follows 

5^ A = eij , S^B = €75^ , 

= ^Y^^{A + ^,B)e + ^{F + ^,G)e, (13) 
SQF = eYd.^, 6'^G = eYlM. 



After eliminating auxiliary fields F and G with the use of their equations of motion the lagrangian 
in Eq. (|T^ becomes 

C = ^d'^Adf^A + ^d^Bd^.B+'-iP^^d^i; 

- ^g'{A' + B^y -zgij{A + ^,B)2p. (14) 

Without auxiliary fields supersymmetric transformations ()14|1 do not form a closed algebra. However, 
introducing light-cone coordinates and defining light-cone fields one can find closed superalgebra for 
"plus" - components of the fields. In light-cone coordinates general four-vector is represented as 
X± = {1/V2){X°±X'^), and X'^Y;, = X+Y_+X_Y^-X''Y\ where i=l,2. For spinors the appropriate 
definition is X± = |7±7=pA, so that A = A+ + A_. Then choosing e = e_ the supersymmetric 
transformations take the following form 

6^ A = e_^+ , 

5^5 = e_75V'+, (15) 
5^^^ = i9+(A-755)7_e_ . 

Now, let us introduce the local conformal Wilson twist-2 operators appearing in this model for 
unpolarized and polarized cases 



(16) 



Ot = iM^d^yi^l.cfC^)^^, (17) 



oi, = A(za,)'+^c;5(|±)A+i?(z9,)'+^c;5(|±)i?, (18) 

= A{^^^r^C]i\{^y + B{^^^r^^^^^^^^ (19) 
where V=d — d,d=d + d and C^{z) - Gegenbauer polynomials 

. (-i)jr(n + .) r(,, + 2.) _ ^ _ 

n\ r(i^) T[2n + 2u) dz^ '- 
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Moreover in this model one can introduce also the following fermionic (by quantum numbers) oper- 
ators 

O'^ = ^+(^5+)'+^P,^?(^)(A + 755), (21) 
= ^+(^5+)'^S5i^^if (^)(A + 75i?), (22) 

where Pi°''^\z) - Jacobi polynomials 

Pt^\z) = ^(1 - z)--{l + z)-P— [(1 - zni + zfH - zY] ■ (23) 

As was already mentioned in introduction, we will be interested in the renormalization properties 
of these operators. It should be noted, that in the singlet case there is mixing between bosonic 
operators formed by fermion and scalar fields. Also, in the non-forward kinematics, in contrast to 
forward case, the operators (fT^ and (fTHj) will mix under renormalization not only with each other, 
but also with the total derivatives of these operators. 

The scale properties of twist-2 conformal operators, we are interested in, could be deduced from 
scale properties of corresponding distribution amplitudes. So, in what follows, we will first determine 
the evolution kernels for DA in the cases of interest. 



2.3 Unpolarized case 
2.3.1 Non-singlet case 

In non-singlet case we need to calculate the evolution kernel of DA, which Gegenbauer moments are 
given by matrix elements of operator sandwiched between quark states. In leading order there 
is only one diagram shown in Fig. ^a and in next-to- leading order there are three diagrams shown 
in Fig. 121 a, Elb and |21c ^. For the diagram evaluation we used the method proposed by Mikhailov 
and Radyushkin P^ ll5j. The calculation details together with diagram by diagram answers may be 
found in Appendices A and B. The final result for the non-singlet non-forward evolution kernel is 
[a 



47r ' 



^Ns{x,y) = £vi.°i(x,i/) + (£)'v«(x,i/) (24) 

Vf)(x,y) = 9{x<y)i-4---\n{x)--\n{x) + -\n(l--] ~-\n(l-- 
{ y y y y \ yj y \ y 



+ -\n\x) + 2-\n(x)\n(y)- -\n\y) + 5(x-y))- + { _ (26) 

y y y } \ y ^ y J 

In the forward case the evolution kernels are given by 

P^^liz) = 2-z-5{l-z) (27) 

p5^).(2) = -Az-2z\n{z) + 2z\Yi{z) + {l + z)\n^{z) + 5{l- z) (28) 

^One need to include also graphs with external self-energies which give ^-function contribution in the coordinate 
space 
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First three terms in the second hne in Eq. (j26j) represent non-symmetric piece of the kernel |/yV/vs(a;, y) 
with respect to the exchange of variables {x ^ y) 

yy^Nsix, y) = xxVNsiv, x) (29) 

and, as was shown in Ref . |14t HE] , lead to non-diagonal part of anomalous dimensions matrix in the 
basis of Gegenbauer polynomials (see below). 



2.3.2 Singlet case 



In singlet case there is mixing between operators (fT^ and (fTH|) and the evolution kernel takes now 
the matrix form 



V(x,?/) 



V^^{x,y) V^^{x,y) 
\i^ix,y) \^ix,y) 

with the elements of this matrix in leading order (diagrams Fig. ^a - Ed) given by 



(30) 



vX^(x,i/) 



vj;^(x,y) 



1- 

'y 





^ X 


'r-^ X\ 


hi 








V y 


) 



9{x <y)<2x 



e{x <y){- 



X ^ X 

y '^y 

X ^ X 

y ^y 



XA^^y) 

and next-to-leading orders (diagrams Figs. 121- El) 

Vij(x,y) 



e{x<y)\-2^Six-y)\ + 



X ^ X 

y *^y 



(31) 



+ 



e{x < y) 

X 



-6- + - In - 

y y \y 



X 



— ln(x) 

y 



X 



y 



X 



-ln(x) + - l-- ln l-- 



In^(x) + 2- ln(x) ln{y) - - In^iy) + 6{x - y) 

y y y 



yi'Jix,y) 



X < y) < —6x — X In ( — ) + 3a; ln(a;) — 3a; ln(a;) — In ( 1 



y 

X ^ X 

y ^y 

X 

y 



X 



(32) 



-|- xln (x) — 2xln(x) ln(?/) + xln (y) 



1-1 


^ X ^ 


X \ 




K y^ 


^y ) 



(33) 



6 1 
- + - 

y y 



X 



y 



y 



+ 



12 1 
ln^(a;) — - ln(x) ln(y) + - In^(y) 



y 



y 



y 



1 


ln(x) 


1 


+ - 


+ — 


y 


yy 


}- 


/ a; ^ 


X \ 




\ y^ 


^y ) 



(34) 



X 



X'^ix.y) = ^^(x<y)<|2 + ln(-) +ln2(a;) + 21n(a;)ln(|/)-ln2(y) + 5(a;-y) 



+ 



X ^ X 

y ^y 



(35) 
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The corresponding evolution kernels in the case of forward scattering have the following form 





P^'\x,y) = { ^ ] (36) 



Pji^W = -7z-2z\n{z) + z\n{z) + {l + z)\n\z) + 6{l-z) (37) 



P^^Jiz) = -7 + z-2\n{z) + {l-2z)\n{z)+\n\z) (38) 
P^^Jiz) = l-7z-2z\n{z) + {z-2)\n{z)-z\n\z) (39) 
P«(^) = z{l-2\n{z)) + 6{l-z) (40) 

2.3.3 Fermionic operator 

The calculation of evolution kernel for fermionic DA is similar to the calculation of evolution kernels 
for bosonic DA. However, in this case there is one subtlety worth to mention. It is necessary to keep 
track of what are the meanings of x and y in the expression for evolution kernels. 

Now, for example if we agreed that x and y are the fractions of momentum for scalar field in 
fermionic operator, then we have the following answers for evolution kernels in leading and next-to- 
leading orders (see diagrams in Fig. ^e and Fig. ^ 



Vf!^r^(x,y) = e{x<y)\-^^\-6{x-y) (41) 



V«(a:,2/) = eix<y)l^2 



y. 

H^) Hy) _ ^ 0. > r in^(^) 
y y S 



1 _ ( 1 f X 

+ - ln(x) - 6{x > y) < - \n ( 1 

y [y \ y 

+ e{x<y)l- + -\n(l--] +-ln(-]] + 5{x-y) (42) 

[y y \ yj y \yj ) 

Evolution kernel in the case of forward scattering has the following form 

= 2-5(1-.) (43) 
pW(z) = -Q-\n{z)-2\n{z)+\n'^{z)+6{l-z) (44) 

2.4 Polarized case 

In polarized case the results for each contributing diagram differ from those in unpolarized case only 
by overall sign. 
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2.4.1 Non-singlet case 

In non-singlet case the diagrams Fig. ^a and Fig. |21b change the overall sign and thus the final result 
for the non-singlet evolution kernel in polarized case has the following form: 

^^Nsix^y) = ^vl,%,y)+[^yYi,%,y) (45) 
^N%,y) = eix<y)U- + -Hx) + -Hx)--(l--)\n(l-- 

[ y y y y \ yj \ y 

+ ^ In^x) + 2^ ln(x) \n{y) - ^ ln'{y) + 5(x - y)| + ^ ^ ^ ^ (47) 

In the forward case the evolution kernels are given by 

Pfsi^) = -2-z-5{l~z) (48) 
Y>^^s{z) = 4:z + 2z\n{z) + 2z\n{z) + {l + z)\n^{z)+5{l-z) (49) 

2.4.2 Singlet case 

Changing the overall sign of diagrams in the cases where it is necessary and taking into account 
mixing between operators (|T7j) and (|T9|) the evolution kernel in the singlet case has the following 
form 

, %^{x,y) %^{x,y) 
V(x,i/) = I ^ ^ I (50) 

\i^{x,y) \A^,y) 

with the elements of this matrix in leading order given by 



^A^h \ nf nT^I f X X 

\i,{x,y) = d[x<y) 



y } \ y ^y 

( 

y ^y 



Vr(^,l/) = 9ix<y)\2-Six-y)] + ( ""^"Z^ ) (51) 
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and in next-to-leading order 

\^,{x,y) = 0(x<y)|6- + -lny +-ln(a:) + -ln(x)--(^l--jln(^l-- 
+ ^ ir^^x) + 2- ln(x) ln(y) - ^In^y) + 5{x - y)\ + (l^l) (^2) 

y y y } \ y ^ y J 

~ (1) f f x\ fx 
Y^^{x,y) = 6{x < y) <.6x + x\n i— j + x\n{x) — xln{x) + In ll 

+ xln2(x)-2xln(x)ln(?/)+xln2(y)|- ^ (53) 

vj)\x,i/) = ^(x < y) |-- - iln ( - ) + iln(x) - -ln(x) - — In ( 1 - - 

[ y y \yj y y yy \ y 

+ iln^(a:)--ln(x)ln(y) + iln^(^)|-f (54) 

y y y ] \ y ^ y J 

%^{x,y) = ^(a;<?/)|-2-ln|^^^ +ln2(x)+21n(x)ln(?/)-ln2(?/) + (5(a;-y) 
/ X ^ X \ 

+ - 55 

\y ^y J 

The corresponding kernels in the case of forward scattering are found to be 



,(0) 



, -2^-5(1-2) -2 

[x,y) = [ . . (56) 

-2z -5{l-z) 



v!pi{z) = 7z + 2z\n{z) + z\ji{z) + {l + z)\n\z) + 6{l-z) (57) 

P^^Jiz) = 7-z + 2\n{z) + {3-2z)\n{z) + \n''{z) (58) 

^l^Jiz) = -l + 7z + 2z\n{z) + {3z-2)\n{z)- zin^iz) (59) 

P^'Jiz) = -z(l + 21n(^))+5(l-z) (60) 

2.4.3 Fermionic operator 

As it can be easily seen from Eqs. (j2H) and (j22j) there is a simple relation between fermionic operators 
in polarized and unpolarized cases 

^IT = -75^17 (61) 

Thus, the evolution kernel for fermionic DA in polarized case differs only by overall sign from unpo- 
larized one. 



3 Confer mal Ward Identity approach 

As was shown in Refs. [T71 ^\ non-diagonal part of anomalous dimension matrix may be obtained 
by studying renormalized Conformal Ward Identities (CWI) for the Green functions of elementary 
fields with the insertions of corresponding conformal operators. To deduce the necessary identities, 
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one starts with the generating functional for renormahzed disconnected Green functions of conformal 
operators. Exploiting its invariance with respect to transformations of collinear conformal group we 

{\0,i]{5''x)) = -{[{^''OM - , (62) 

where < A > denotes averaging over the vacuum of the time ordered product TAexp(i[S']) and 
X = Hi 0j stands for the product of elementary fields at different space-time points. The left hand 
side of the above equality is the differential operator acting on renormalized Green function {[Oji\x))-i 
which thus should be finite. Hence, the right hand side of the equality is also finite. Using this fact 
we may determine the scale and special conformal anomalous matrices for conformal operators under 
consideration. 

Now, the crucial point is to note, that scale and special conformal anomalies are in fact related. 
The relation between matrices of scale and special conformal anomalies could be obtained by consid- 
ering the action of the commutator of special conformal transformation with dilation [5^ ^ 5^] = 6^ 
on Green function with conformal operator insertion. As a result we obtain the following relation 



a(/)+r(«,0 + ^K0,7(«) 
a 



0, (63) 



where ^'^{a, I) is matrix of special conformal anomalies, a(/), b{l) is some matrix coefficients and 7(a) 
are the anomalous dimensions matrix, we want to find. Due to the fact, that matrices a{l) and b{l) 
do not depend on coupling constant (in contrast with all other quantities, appearing in Eq. (jHSl)), 
it is possible to solve this equation recursively at a given order of perturbation theory and find the 
relation between non-diagonal elements of anomalous dimension matrix 7^^^ in n'^^-order and matrix 
of special conformal anomalies 7^(a, /) in (n — l)-order of perturbation theory. For example, at LO 
we have 7^° = Sjk'Jk- It should be noted, that in the absence of gauge field the final expression 
is considerably simplified and one only needs to know the matrices a{l) b{l). More details about 
this method in the case of fermionic operator could be found in Appendix C. Below we present the 
results, obtained within this framework, for non-diagonal parts of anomalous dimensions matrix of 
conformal operators, we are interested in. 

3.1 Unpolarized case 
3.1.1 Non-singlet case 

Consider first the case of non-singlet bosonic operator. Solving Eq. ()63p in the next-to-leading order, 
it is easy to find the following general expression for non-diagonal anomalous dimensions matrix 

[13 CHI 



-^^r)rf„.(6o + 7*-^^r) , (64) 
where 



In 



n+l n+2 

is the corresponding leading order anomalous dimension of operator p6|l in forward case and 

, + ' „ if ,-*>0 and even 

otherwise 
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Taking into account the next-to-leading order anomalous dimension in forward case 

2 2 4 2 2 4 



In 



n + 1] 



(n + 1) 



+ 



n 



'n + 2y 



+ 



+ 



[n 



2)2 n + 2 



2 Si (n+1) 2SAn + 2) ^ 
+ ^ + 1 



n+l n+2 
it not hard to verify that the relation (jHJ holds true. 



(66) 



3.1.2 Singlet case 

In the singlet case the expression for non-diagonal anomalous dimensions matrix should be modified to 
take into account operator mixing. In the case of Wess-Zumino model the corresponding expressions 
are 



ij,k 
Mm 

Mm 

</,0 ND 



^^if - rf,,, (bo - ) - {^^if - ) 



(67) 
(68) 
(69) 
(70) 



where "'^'jf'^ are leading order anomalous dimensions of singlet operators p6|) and 1)18^ : 



/n /n 



Wo/ 
/ra /n 



2 n+1 



n+1 n+2 

2 2 
n + 2n + 1 



n 



1 2 



(71) 



Note, that in our normalization of Gegenbauer polynomials, the anomalous dimensions of conformal 
operators (fTI)|) and (fT^ in singlet case differ from corresponding anomalous dimensions of the usual 
(non-conformal) Wilson twist-2 operators appearing in the description of DIS by shift in argument 
n n + 1 and normalization factors. Here and in what follows we will show explicitly additional 
normalization factors or For example, our results for ^^•jn and '^^'jn in Eq- f)71|) contain 
explicit multiplication by these factors. 

The next-to-leading order diagonal elements of anomalous dimensions matrix in this case are 
given by 



In 



In 



In 



In 



1 



(n+1] 



+ 



1 



n 



n 



'n + 2) 



(n + 2) 



+ 



7 



n + 2 



+ 



2 Si (n + 1) 2Si(n + 2) 



n + 1 
2 



n + 2 
1 



(n + l)3 (n + l)2 
2 1 



+ 1 
7 

n + 1 
2 



+ 



[n + 2)2 
1 



1 2Si(n + l)\n + l 
+ ^ + ^ 



n + 2 



n + 1 



7 ^2Si(n + l)\ 2 



n 



1)2 n + 1 (n + 2)3 (n + 2)2 n + 2 n + 2 J n + 1 



(n + 1) 



+ 



n+1 (n + 2)2 n + 2 



+ 1 



(72) 
(73) 
(74) 
(75) 
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In singlet case the relation similar to Eq. (jS)) has the following matrix form 







dx ( C'J\x-x) Cy+i(x-x) ) 



n (^^ln,k ^^ln,k\ ( Cl'\y - y) \ 

= E (76) 

It could be easily verified substituting explicit expressions for NLO evolution kernels Eqs. ()32H35p 
and NLO anomalous dimensions Eqs. ()72II75|) . 

3.1.3 Fermionic operator 

The expression for non-diagonal anomalous dimensions matrix of operators with fermionic quantum 
numbers (pT|l could be derived similar to the case of bosonic operators (for more details see Appendix 
C). It is given by 

%T = (N°^-'^^7r)4:^^(/5o-^-7r) , (77) 

where 

fer^(0) = (_l)n+l 2 _^ (78) 

n+1 

is the corresponding leading order anomalous dimension of fermionic operator (j2ip in forward case. 
The matrix is defined through the derivative of corresponding Jacobi polynomial over its indices: 



d 

dv 

jfer 



3 

jfer 0(1-0), 



= 2E4l^i W (79) 

^^=0 fc=0 



The explicit expression for d^J. {j > k) has the following form 

= l(l±i)±(zi)!^(^^ (80) 

^'"^^j-A; J + k + 2 J + 1 ^ ^ 

Taking into account the next-to-leading order anomalous dimension of this operator in forward case 

'^^ {n + iy^ ^ ' \{n + iy n + 1 n + 1 J ^ ' 

it easy to verify, that the relation similar to Eq. (jH)) 

f'dxP^'''\x - x)VUx,y\a) = Y,''''lAa)Pt'^'\y - y) (82) 

"^0 fc=0 

is valid. 



3.2 Polarized case 

In the polarized case non-diagonal anomalous dimensions matrices have the same expression as in 
unpolarized case (see. Eqs. (jMj) . (pTjl - lfTOj) and (fTTjl ). The only change, which should be done, is the 
replacement of unpolarized forward anomalous dimensions 7j with polarized ones 7j. 
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3.2.1 Non-singlet case 



Taking into account the change of sign in diagrams in Fig. ^a and Fig. |21b we find that polarized 
leading order anomalous dimension is given by 



In 



whereas in next-to-leading order we have 

2 2 



+ 



n+ 1 n+2 



1 



(83) 



In 



+ 



+ 



+ 



{n + lf {n + lf n + l {n + 2f {n + 2f n + 2 
2Si{n + l) 2Si{n + 2) 



n+l n+2 
It easy to see , that the relation (jH)) holds true. 



+ 1 



^4) 



3.2.2 Singlet case 

In singlet case the leading order anomalous dimensions are 

2 2 



wz 



In 


In 


in 


In 



n+l n+2 

2 2 

n+2n+l 



2 n+l 
n+l 2 

-1 



^5) 



In the next-to leading order the diagonal elements of the anomalous dimensions matrix have the 
following expressions 



in 



In 



In 



In 



+ 



+ 



+ 



7 



+ (n + l)2 n + l (n + 2)3 {n + 2y n + 2 

2^1 (n + l) 2Si{n + 2) 



n + l 
2 

[n + l)^ 
2 



n + 2 
3 7 



+ 1 



+ 



2^i(n + 1)\ n + l 



n + l)2 n + l (n + 2y n + 2 



n + l 



+ 



7 2Si{n + l] 



(n + l)2 n + l (n + 2)3 (n + 2)2 n + 2 

2 12 1 

+ ^ + 1 



n + 2 



2 
2 



n+l 



(n + l)2 n + l (n + 2)2 n + 2 

4 Supersymmetric Ward Identity 



(86) 
(87) 
(88) 
(89) 



Now, let us recall, that all operators considered above form the representation of = 1 supersymme- 
try algebra. Similar to CWI we may consider the supersymmetric Ward identity [211 (i^^ derivation 
goes along the same lines as the derivation of CWI). As we will see in a moment, there is a host of 
remarkable consequences for anomalous dimensions of our operators following from supersymmetric 
Ward identity. To begin with, let us introduce the following combinations of conformal operators for 
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unpolarized (Eqs. (HHD, (UHl) and (jUD) and polarized (Eqs. (HZD, (UHl) and (j22I)) operators. 



s], 

V,.i 



1^ I q4 
2 



J + 1 



J + 2 J + 

fer 



(90) 

(91) 
(92) 



where the coefficients in front of operators define the matrix, diagonahzing the matrices of forward 
anomalous dimensions (f7T|) and (jH3j) . 

Under restricted light-cone supersymmetry transformations the combinations of conformal oper- 



ators introduced above transform as follows 






s^'sh = (1 - (-i)^>-vj_,, 




(93) 


s'^si = {i-{-iy)evl 




(94) 


S^Vl = {1 + {-iy)eVl,, 




(95) 






(96) 






(97) 


The transformation low for the operators Uj^k foolowing 


from the observation Uj^k = 


— TsVj^fc and the 



Eq. (jnni)- Now, we see that it is natural to combine bosonic and fermionic conformal operators into 
= 1 chiral supermultiplet. 



(9^ 



where Vj-i^i is the Majorana fermion build from the Weyl spinor X . 

It is easy to find, that the renormalized supersymmetric Ward identity in the regularization 
scheme, preserving supersymmetry, has the following form {Sj^i denotes vector of operators 5"^^ 



and Sli) 



{[S,,]5'^X) = -{5'^[S,i]X) - {i[S,i]{5'^S)X) and {5^[S,i]X) = finite 



(99) 



where we used the fact, that renormalized action in supersymmetric regularization is invariant with 
respect to supersymmetry transformations {i[S ji\{5'^ S) X) = 0. The operators and mix 
under renormalization and thus we define renormalized operators as (square brackets correspond to 
renormalized quantities) 



E 

A;=0 



lly 12y 



jk 











(100) 



kl 



and as a consequence the renormalization group equation for these operators is given by 



d 



(iln/i 



51 



31 



E 

fe=0 



2y 22^5 



jk 



51 
52 



(101) 



kl 



Now, from supersymmetric Ward identity we get (cTfc = |(1 — (—1)'^) and Zjk = for /c > j) 

3 k 



EE 

fc=0 fc'=0 



lly 12y 
21y 22y 



CTk 



jk 



{Zx>^}k-1,^ 

{Z^ }kk' 



k'l\ 



finite 



(102) 



15 



1/e poles in (|1U2|) cancel, provided 



t{ 

k=0 
j 

E 

fe=0 



jk 



jk 



Ok [V> 



k-lA. 



Ok [Vk- 



1,1 



fc=0 

j 



jk 



Ok m 



kl\ 



k=0 



jk 



Ok [Vl 



kl\ 



a. 



a. 



E 

fc=0 

j 

E 

A;=0 



[V; 



kll 



jk 



kll 



(103) 
(104) 



Noting, that Uj^k = — TsVj^fc, it is easy to derive analogous relations for polarized operators Vj^k- 
Taking into account linear independence of fermionic operators [Vki] we finally get the following 
relations ^T] 



/2n+l,2m+l 

I2n+l,2m+l 

2V 

/2n+l,2m+l 

22 "5 

72n+l,2m+l 



- 22 

/2?i,2m 

- 2U^ 

/2n,2m+2 

- 12^^ 

/2n+2,2m 

~ /2n+2,2m+2 



fer 



fer 



fer. 



fer 



T2?i,2m; 
72n,2m+l5 
'72n+l,2m, 
72n+l,2m+l, 



m < n, 
m < n — 1, 
m < n , 
m < n , 



and 



12 5 

72n+l,2n+l 







12^^ 
/2n,2n 



. 



(105) 
(106) 
(107) 
(108) 

(109) 



These relations allow us to find anomalous dimensions of bosonic operators from known anomalous 
dimensions of fermionic operators and vice versa. In particular, substituting explicit expression for 
non-diagonal anomalous dimension of bosonic operators, one can determine matrix dj^k for fermionic 
operators Eq. ()80j) . The found relations between anomalous dimensions of conformal operators is 
special case of existing super conformal relations for anomalous dimensions of conformal operators in 
super symmetric theories. 

To obtain from the results above the corresponding relations for anomalous dimensions of confor- 
mal operators considered earlier, one will need the following transition formulae among anomalous 
dimensions of conformal operators (jnn|) and (PT|) and anomalous dimensions of operators (fTI)|) and 



m 



( 



k+\ 
1 

/c + 2 
1 

k+\ 
1 

V k + 2 



7 



12, 



7 



21 



7 



22, 



7 



\ 




/ fc + 2 


k + 2 


2 






J + 1 
k + \ 


2 

k + \ 


J + 1 
2 




1 


J + 1 


2 


J + 1 




2A; + 3 


k + 2 


A; + 2 


2 






J + 2 
k + \ 


2 

k + \ 


J + 2 
2 


' jk 




\ J + 2 


2 


J + 2 



1 



\ j 



(110) 



jk 



Now, let us explore some particular limits of equations ()105|) - (ll09p . From Eq. ()109|) it follows, that 
(here and below 7^^ = 7^, 7 = 7^, 7 = 7"^ and indexes S {V) correspond to unpolarized (polarized) 
cases) : 

7-1-1 2 



^^7i + 



7, 



;iii) 



2 J + 1 

where j is odd for unpolarized case and even for polarized case. This equation is analogous (up to 
redefinition of anomalous dimensions) to the well-known Dokshitzer relation (last paper in Ref. [2]), 
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which was found first empirically in LO in the case of supersymmetric QCD (the anomalous dimen- 
sions in supersymmetric QCD could be obtained from those in ordinary QCD through substitutions 
Ca = iV„ Cf = and Nj = \N,.). 

From p05|) and p08|) . we also have (for odd j) 



7,- 



J • 2 



7 



+ 



J + 2 



7j+i 



2, 

j 
2 



'7 — 1 

3 ^ 



fer. 



7j-i 



fer 



J + 1 



7j 



(112) 
(113) 



Moreover, from equations ()106|) and pOTj) we get the following relations between diagonal and non- 
diagonal elements of anomalous dimensions matrices of conformal operators 



7,- 



J + 1 



0V 



7,- 



J + 1 



la'' 



fer. 



-1 



3 + 2 2 '-^ 2j + 

where j is odd for i = S and k = V, j is even for i = V and k = S and 

_ J i/ii/) . « , J -^»'' « 2 



7ij-i 



A 



J+2 



J + 2 



7i+i,i-i 7j+i,j-i- 



(114) 



(115) 



In the leading order A*^_^j_^ is equal to zero, but starting from next-to-leading order it acquires 
nonzero contribution. Here, we would like to note, that all these relations are valid in all orders of 
perturbation theory provided anomalous dimensions were evaluated in supersymmetric scheme. In 
Wess-Zumino model such scheme could be identified with the usual MS-scheme. 

Next, the forward anomalous dimensions matrix in singlet case ()7ip is diagonalized with the help 
of matrix D 



D 



so, that in leading order we have 

[DrwzD-i] 



11(0) 



J + 1 
2 

"jT2 



J + 1 




(116) 



J + 2 



(117) 



In next-to-leading order of perturbation theory the elements of forward anomalous dimensions matrix 
of bosonic operators are given by (j72p - (j75|) . Using the same diagonalizing matrix D (jll6|) . as in 
leading order, we get 



[DFwzD-i] 



11(1) 



:ii8) 



(J 



+ 



1 



6 , 2 5i(j + l) 



(j + If 
4 



J + 1 
4 



J + 1 



+ 1 



1 



J+1 J+2 



+ 



6 2 5i(j + 2 



(j + 2)3 (j + 2)2 J + 2 



J + 2 
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So, in leading order the following relations between elements of anomalous dimensions matrix are 
valid 



^ 7i + -^^^Ij 



J + 1 
2 



7i 



7,+ 



J + 1 J + 1 



7i 

7,+ 



7, 



(119) 
(120) 



2 J + 2 

However, in next-to-leading order the relation p20|) does not hold, while the relation pi9p is valid 
to all orders of perturbation theory. Note, also, that from the relation ()119p it follows that 

2 ,,,, j + 1, 



7j 



7i 



7, + ""7i 



(121) 



J + 1 " 2 

which is also valid to all orders in perturbation theory and in leading order coincides with the 
relation (|12(jp . Moreover, modifying the diagonalizing matrix (|116p it is possible to show, that the 
forward anomalous dimensions matrix will remain diagonal also in next-to-leading order, and, what 
is more important, in all orders of perturbation theory. Such matrix D is given by 

2 



2. 



D 



J + 1 

7j 



1 



■4,4,, 



(122) 



In polarized case the forward anomalous dimensions matrix (jHK|l may be diagonalized with the help 
of the same matrix D ()116|) . as in unpolarized case and it is possible to find, that 

2 



DFwzD 



-1 



(0) 



J + 1 




1 







J + 2 



1. 



(123) 



In next-to-leading order the elements of forward anomalous dimensions matrix in polarized case are 
given by Eqs. (jHUjl - lfH^ . Using the same diagonalization matrix, as in leading order, we get 



DFwzD 



-1 



(1) 



(124) 



6 , 2^i(j + l) _^ 



(j + l)3 (j + l)2 J + 1 J + 1 

4 4 8 2 1 6 2 5i(j + 2) 

+ — - 1 



J + 1 J + 2 (j + 2)2 (j+2)3 (j + 2)2 J + 2 J + 2 

At LO we have three independent relations among elements of forward anomalous dimensions matrix: 
Eq. (fTTT]) together with Eq. (fTT^ . Eq. (|TT^ and Eq. where at LO we have ^'j+ij-i = 0. 

Therefore, we may conclude, that in forward case in leading order of perturbation theory, it is 
sufficient to know only one diagonal element of anomalous dimensions matrix (|7T|) or anomalous 
dimension of fermionic operator ^^'^Jj to reconstruct the whole matrix. 

At NLO the righthand side of Eq. (jll4|) is nonzero and as a consequence matrices (jll8j) and (jl24j) 
become triangular. So, to reconstruct the forward anomalous dimensions matrix one needs to know 
additionally any other element of forward anomalous dimensions matrix or non-diagonal anomalous 
dimension of fermionic operator *^^'"7j 

^We thank Anatoly Kotikov for fruitful discussions of this issue. 
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5 Solution of NLO evolution equation 



Our primary goal is the construction of multiplicatively renormalized conformal operators up to 
next-to-leading order in perturbation theory. The non-singlet case is simple and the corresponding 
solution could be easily found in full analogy with 0^-theory 14j. So, here we will be interested 
only in singlet case with bosonic operator mixing. In previous section we shown, that the twist-2 
conformal operators form a chiral supermultiplet of = 1 supersymmetry algebra and derived their 
transformation rules under supersymmetry. It turns out, that in the singlet case, it is sufficient to find 
only the expression for multiplicatively renormalized fermionic operator. The multiplicatively renor- 
malized bosonic operators up to next-to-leading order in perturbation theory could be obtained then 
via supersymmetry transformations from the expression for multiplicatively renormalized fermionic 
operator. 

As we already mentioned before there is a close relation between DA and conformal operators, 
we are studying. To see it explicitly, let us consider the case of non-singlet bosonic operator. The 
eigenvalues of the anomalous dimensions matrix 7 of this operator are 7^ = jjj. Therefore, the 
renormalization group analysis provides for multiplicatively renormalized operators Oji: 

d,i{i?) = exp (^j'^ jj, {g{t))^ OM)- (125) 

Moreover, the operators Oji could be completely expressed through Oji\ 

k 

Ojiifi') = Y,B,,{g{f,'))dM{^'). (126) 

k=0 

Then, the evolution of distribution amplitude is given by: 

<P{x,Q') = ^0„(a;,a(Q^)) exp / -7n(^?(t)) (0|O„„(Qg)|P)-^ (127) 



where 



n=0 V '^O 



00 _ 
XX 



[x, a{Q')) = J2 -ufl"^"^ - ^)BknHQ')) ' (128) 



k=n 

where A4 is some normalization factor (see below). So, in this section we will first determine the 
eigenfunctions of the fermionic evolution kernel. Then, it is straightforward to write down the mul- 
tiplicatively renormalized fermionic operators. Our treatment here is close to the work of Mikhailov 
and Radyushkin on 0^-theory in 6-dimensional space-time [TH I16j. 

In leading order of perturbation theory the evolution equation Q for the fermionic operator ()21|] 
could be easily solved employing separation of variables. 



[x,n 



'1lo = E«" l^X^ xP(i'°)(x-x), (129) 

n=0 ^ ^ 



where ^^'^'jn'^ is the LO anomalous dimension of the fermionic operator (j2H) and xPji'''^\x — x) = ipn{x) 
is the solution of the eigenvalue equation 

vo{x, y)'4)n{y)dy = A„^/'„(x) , (130) 
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where A„ = (— + 1) and Vq is defined as 

Vfi°^ {x, y) = v,{x, y) + 6{x - y) (131) 

The easiest way to find ipnS is to note, that the combination yvQ{x,y) = v{x,y) is symmetric under 
the change x y: v{x,y) = v{y,x). Another property of Vo{x,y) is that its convolution with the 
polynomial of degree is the polynomial of degree M with M < N. Thus, the functions ilJn{x)/x 
should be the polynomials orthogonal to each other on [0, 1] with measure x. 

Now, let us proceed with NLO evolution kernel. It is convenient to rewrite the resulting expression 
Eq. ()42|) for NLO evolution kernel of the fermionic operators in the following form 

V}^^ = eix > y) [-\nix) - - Inf 1 - - ) 1 + 0(x < y) [-\n{y - x) + - 

yy y \ yj) U y 



\n\y) 2 ln(x) ln(?/) ) [ In 



+ e{x<y){ ^ + ^"^ |+g(x>y)|— (132) 

+ e{x < y) + 5{x - y) 



It is easy to verify, that the combination y^^^} is symmetric with respect to transformation x 



y if one takes into account only the first line of Eq. ()132|) . The second and third lines of the 
expression above constitute non-symmetric piece and thus lead to the non-diagonal terms in the 
Jacobi polynomial basis (excluding 5-function, of course). To find a solution of NLO evolution 
equation we will need the expression of NLO evolution kernel in the basis of LO eigenf unctions. To 
simplify further notation it is convenient to introduce following definitions: 

= "(i+/--H^-^) (_i)»p(i..^)(,_^) (133) 

where J\fn{a,(3) = 2°+/^+^ is the normalization factor and w{a,(3,x) = 

(1 — x)"(l + x)^ is weight function. These are just the solutions of the eigenvalue equation ()130p and 
its generalization: 

/ dyV^'\x,y)^';}f{y) = ^ ^'^fiy) , (135) 

Jo n -\- I -\- u 

where 

Vj'\x,y) = e{x<y)-(-) (136) 

y \yj 

To write down the evolution kernel in the basis of LO eigenf unctions, it is convenient to express the 
NLO kernel as the convolution of LO evolution kernels and their derivatives: 

vo{x,y) = d{x<y)^- (137) 



Vo{x,y) = (AvW(x,y) 



v^,.. , e{x<y)-\n{^] (138) 

Then it is an easy exercise to derive, that: 

V'^^\x, y) = -2vo{x, z) ® vo{z, y) + Qvo{x, y) + 2vq{x, y) + m(x, y) , (139) 
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where ® denotes the convolution and 



X 



u(x, y) = 6(x > y) < — \n(x) — - In 1 > + 6(x < y) < — \n(y — x 

[y y \ yj) [y 

It is not hard to show that u{x,y) is diagonal in the basis of LO eigenfunctions: 



(140) 



dyu{x,y)¥^''\y) = {-ir^' 



(n + l) 



1 ^2S^{n + l] 



n + l 



(141) 



Now, what is left is to specify the action of dot kernels on LO eigenfunctions. As a starting point we 
take the following eigenvalue equation 



Vj'\x,y)(^¥jf{y) 



-1) 



71+1 



X) 



n + 1 + u 

Differentiating this equation with respect to u and putting it afterwards to zero we get 

d 



(142) 



Vo{x,y) ® ^^^'^\y)+Vo{x,y) ® — (^i^f (^/)) 



u=0 



(143) 



So, to find an action of dot kernels on LO eigenfunctions we need to know an expression for the u 
derivative of The latter is given by 

d 



= 2 dkn'^k^'^^ (x) + diagonal part , 

^=0 k>n 



(144) 



where 



d 



■jk 



1 {j + i) + {-iy-\k + i)k + i 



j -k j + k + 2 



J + 1 



(145) 



We failed to find diagonal part in the equation above, however it is quite easy to derive diagonal part 
directly for Vf){x,y) ® \l/n'°^(?/). In this way, we get from Eq. ()143j) 

(_-\\n+l °° 

M^,y) ® ^t'\y) = u;iTY^n''\^) -J^^Uif - (i46) 



k>n 



We can summarized the actions of building blocks of NLO evolution kernel on LO eigenfunctions as 
follows 



Vo{x,y)®<i/il^^\y) 
vo{x,y)(E)^^n'''\y) 
u{x,y)®¥^^^\y) 



n+l " ^ ' 
-1 



in+l 



k>n 



\n+l 



1 ^2SU^| 

(n + 1)^ n + l J 



(147) 
(148) 

(149) 
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With all these formula at hand it is easy to derive that 

-2{;o(x,z)®t;o(z,y)®vl/a,o)(^) ^ M^fi){x) (150) 



0<J 



k>n 

Note, that substituting the expression for NLO evolution kernel in the basis of LO eigenfunction into 
Eq. (j82p we can easily determine analytical expression for non-diagonal part of anomalous dimensions 
matrix of fermionic operator. 

Now we are ready to write down the solution of ER-BL evolution equation for fermionic operator 
()21|) in next-to-leading order. Note, first, that if the next-to-leading order kernel V^^\x^ y) is diagonal 
in the Jacobi polynomials basis, then the solution of evolution equation Q would be 

=a„(/xg)exp|^-^' ^-^^^{g(t))^ xPi'^'\x - x), (151) 

where *'^''7„(^) = (^) ^""^^i^n^ + [f^f^'^^i^n^ with ^""^-^^n^ and ^""^-f^n^ given by Eqs. and (jHH) corre- 
spondingly. In the case when V'^^\x^ y) contains the non-diagonal terms, the solution of the evolution 
equation will differ from by 0{a) terms. Therefore, we will look for the solution in the following 

form 

^n=[l + ^w)®<pf^^\ (152) 

Consider first a more simple case 6o = 0. Substituting p52|) into the evolution equation and using 
the explicit form of the total kernel ()139p one can find, that in this case W = W{x, y) should satisfy 
the following equation 

[vo,W]_+vo(S)i-f^-vo) = 0. (153) 

This equation could be easily solved noticing, that it is similar to eigenvalue equation for the vi°^ 
kernel written up to O(z^) terms 

[vo,uj]_ + Vo = 0, (154) 
where lj stands for the generator of i/-shifts: 



^(1,0) ^ ^(1,0) , ^A^(i.o) 

u,n n ' u,n 



+ C>(z/2) = {l + uuj)® ^i^'") + 0(z/2) (155) 

u=0 



As a result, we have the following expression for W: 

W = ~co0ivo- 7^) (156) 

and as a consequence we have 

0n = (l + £w^)®^^°)=(l-£(A„-7^)-)®v[;a.o) 

= ^i!;°H^)L_^,(o)(l + 0K)) (157) 

47r V / 

Then, the multiplicatively renormalized operators in this {bo = 0) case could be written through the 
coefficients dkn fll45j) of the expansion ()144|1 : 

On = On + ^Yl dknlTOn (158) 
n>k 
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In the &o 7^ case the equation for the rotation kernel W has even more comphcated structure than 

Eq. (USHD 

boW + [vo, W]_ + vo® {bo + 7^ - ^^o) = 0. (159) 

The effective method to solve this equation is to use matrix representation for the kernels in the 
Jacobi polynomial basis. It then follows, that W ® "^n'^^ has the following form: 

W ® = {bo + 7$ - A„) (A„ -Vo- bo)-' ^vo^ (160) 

Now using the expression for the dot kernel in the Jacobi polynomial basis derived previously, we 
get 

W ^ ^^°) = {bo + -7f ) (Pn^^:^'^ + bo± (0) t (0) , ^""^l ' (161) 

V k>n Tk - In ' -bo 



where Pn is the projection operator Pn'^n'^^ = (1 — 5„fe)^n subtracting the diagonal part. The 
multiplicatively renormalized operators in this case are given by 

On = On + -{bo + ^^^7i°^) ^ J" ^ O,. (162) 



n>k 



6 Conclusion 

It follows from our analysis, that to find multiplicatively renormalized twist-2 conformal operators 
in Wess-Zumino model it is sufficient to find multiplicatively renormalized operators only for one 
member of operator supermultiplet. It is convenient in this case to find a solution for fermionic op- 
erator. For the latter the problem of finding multiplicatively renormalized operators is much simpler 
compared to singlet bosonic operators, where we have operator mixing. All other multiplicatively 
renormalized operators could be then obtained from fermionic operator with the use of supersymme- 
try transformations. Moreover, we found, that the knowledge of fermionic diagonal and non-diagonal 
anomalous dimensions matrices allows us completely reconstruct the forward anomalous dimensions 
matrix in singlet case. Also, the analysis was performed in a simplified supersymmetric Wess-Zumino 
model, we think that most of our findings could be straightforwardly applied to field theories with 
more supersymmetries. 
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Appendix A 

Here we give the diagram by diagram results. The whole calculation was automized with the help of 
Feynman diagram analyzer DIANA [21], calling QGRAF [21] for the Feynman diagram generation, 
and computer algebra system FORM 

The results for each diagram include contributions of this diagram to forward anomalous dimen- 
sion 7m, to forward evolution kernel P{z) and to non-forward evolution kernel \{x,y). Here and 
below, the contributions of diagrams marked by star change the overall sign in polarized case. 
The tadpole diagrams (similar to^d) give zero contribution in forward case and are proportional to 
constant in non-forward case. So, in what follows we do not write explicit expressions for contribu- 
tions of these diagrams. Also, graphs containing external self-energies, giving 5-function contribution 
to evolution kernels, were taken into account in calculation, but not shown explicitly here. 
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Figure 1: One- loop diagrams 
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(a) (b) (c) (d) 

Figure 2: Two- loop diagrams for quark operator Eq. sandwiched between quark states 
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Figure 3: Two-loop diagrams for quark operator Eq. sandwiched between scalar states 
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(a) (b) (c) (d) 

Figure 4: Two-loop diagrams for scalar operator Eq. (|18|) sandwiched between quark states 
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Figure 5: Two-loop diagrams for scalar operator Eq. (|18|) sandwiched between scalar states 
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Figure 6: Two-loop diagrams for operator Eq. ()2ip with fermion quantum numbers 
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Appendix B 



Here we present details on the evaluation of momentum integrals present in the calculation of evo- 
lution kernels of ER-BL equations. Our consideration follows the original work of Mikhailov and 
Radyushkin [TH] ^ 

The main complication of the integrals involved compared to usual propagator-type integrals (for 
which a lot of powerful technics were developed over the past years) is the presence of (5-function 
related to effective vertex. To derive all the formula we used a-presentation for momentum integrals 
as was originally proposed in pH]- For propagators and (5-function insertion we have: 



J daa^ ^ ■ expl^a^fc^ + I (163) 



\ X — - — = — / aaexp <ta\x — - — > (lo4j 



k ■ n\ 1 

In what follows we used the same integration measure convention as used by ^Hl, that is we use MS 
prescription. Technically, it could be accomplished by multiplying each integration d'^k by additional 
factor M{e) = (47r)~'^r(l — e). At one-loop level all integrals obtained after projecting out Dirac 
structure of the conformal operator under consideration and expanding scalar products of momenta 
in the numerator over the denominator factors could be reduced to the following two most general 
integrals 

{27tY [{k - aPy - mf]"i [{k - hPf - mi]"^ 



(27r)'^r(ni)r(n2) 



T(ni+n2--^ / t;aa"i-i(l -a)"2-^5(x-aa-fe(l -a))A5-"i-"2 

(165) 



and 



{27iY [P - m?]"i[(A; - aPy - mlY''^[{k - hPf - m|]"3 



(2.)^r(n0r(n.)r(n/("^ + "^ + "^-2; 

X / {d/32d/33}(l-/32-/33)"^-'/32'"'/33'"''^(a;-«/32-&/93)i?^^"^-"^-"% (166) 

where the following notation is introduced: c? = 4 — 2^, {dii^dii??^ = 0(1 — (32 — P3)dj32dj33 and 

A = -amj-{l-a)ml (167) 
B = -I32ml-(3sml-{l-l32-(33)ml (168) 

The masses ml, m\, m\ provide infrared regularization in the cases where it is necessary. 

On the other hand these expressions for one-loop integrals could be used as building blocks to 
obtain expressions for necessary two-loop integrals. Here we would like to note that evolution kernels, 



•^It should be noted that expressions for J{z) and S{z) functions (to be introduced below) have misprints in [T^ . 
We thank Sergei Mikhailov for fruitful discussions on this subject and performed crosschecks. 
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in which we are interested, are being obtained as a result of apphcation i^i?'-operation to diagrams 
describing renormahzation of effective vertices corresponding to multiparticle distribution functions. 
So, in what follows we give answers for all necessary two-loop integrals before and after application 
of i('i?'-operation. Note, that it would be incorrect to use just i^'i^'-subtracted scalar integrals as 
extra e terms could follow from Dirac algebra for 7-matrices and procedure used for subtracting 
sub divergences provided by /^'-operation should take them into account. 

To perform the calculation described in this paper the following general formulas for two-loop 
scalar integrals are sufficient 



V{h\c,d\l) 



-(4 



7r 



1 



(/ - cP)2(/ - dPy J (k- l f{k - hPy (27r; 



5{x-^) d'^ld'^k 

2d 



(169) 



Before i^Ti^'-operation we have 
V{h\c,d\l) ^ 



-xi)d{xi > yi)J 



+ e{xi)6{xi < yi + y2)J 
Here we have introduced the notation 



Xi 

yi 

Xi 



yi + y2 



-xi)9{xi > yi + y2)J 
0{xi)e{xi < yi)J 



Xi 



yi + y2 



(170) 



d'^k 

J{z) 
J{z) 



M{e)d'^k, xi = X — b, yi = d — b, y2 = c — d > 

11 ~ 

— ln2;H — ln^2;, S(z) 

e 2 ' ^ ' 



II- z I- z 1 + 2, 
2 Inz 



z z 
I- z 1 + z 



\n.z 



e z z z 

The result for \^(6|c, (i|l) after i^'i^'-operation could be obtained from the result without KR'- 
operation performed with the substitution J{z) J{z). 



V b 



c, d 



I ■ n 



Before i^'i?'-operation we have 



-(4 



I ■ n/P ■ n 



(/ - cPf{l - dPf J {k- iy{k - bPy {2tt] 



H^-T^) dWk 

2d 



(171) 



V b 



c, d 



I ■ n 



P-n 



1 



[xi > yi)d{-xi) 



2^2/2 

- 9{xi > yi + y2)d{-xi] 



Xib — 

Vi 



xiS 



Xi 



+ 9{xi < yi + y2)0{xi 



- 9{xi < yi)9{xi] 



xiS 



xiS 

Xi 

yi 



yi + 2/2 

Xi 



bJ{^-l 

yi 
bj 



Xi 



yi + y2 



+ bJ 



Vi + 2/2 

Xi 



yi + 2/2 



+ bJ 



(172) 
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The result for V (b\c,d\-p^^ after i^i?'-operation could be obtained from the result without KR'- 
operation performed with the substitutions: J{z) J{z) and S{z) S{z) 




W {a,b 



(/ ■ n/P-n) 



-(4vr) 



X 



{k - apy{k - hpy 
1 r 2 



{l-kf{l-cPf \ {I -n/P-n) j (27r) 



d'^kd'^l 

\2d 



Before i^'i?'-operation we have 



Wla.h 



(/ ■ n/P - n) 

2e \ \{xi + y + l) + h j yi+y2 

Xi 



1 

- + 2 

e 



0{xi)e{yi + y2 > xi) 

+ 9ix^)e{y^ > x^) - 9{x,)9iy, > x,)e{x^ > y, + y^) In ^'j ~ 

+ e{x, > yi)e{y, + y2> x,) In [^^ ^^j ~ - ^^(-0:1)^(0:1 > y^) In ^ 

[1/1 + I/2J - 1/1 



a-b^c ]_ 

2e\ lixi + yi + b) 



9{-xi)9{xi > yi) + 9{xi)9{yi > xi] 



1 1 



(173) 



(174) 



where xi = x — b, yi = c — b, y2 = a — c, yi + y2 = a — b > 0, x > 0. After .ft'.R'-operation was 
performed the integral above takes the following form 



W {a,b 

a ^ b 



1 

(/ ■ n/P - n) 

1 r 1 



2e 1 \{xi+y + l) + b } yi+y2 



9{xi)9{yi + y2 > xi) 



+ 2 



+ 9{x,)9{y, > X,) In ^ - 9{x^)9{y, > x^)9{xr > y^ + 1/2) In [^^ | ^^j 

yi [yi + 2/2J - yi 



+ 9{xi > yi)9{yi + 1/2 > Xi) In 
a — b^c J_ 



\y\ + 2/2] - xx 



-xx)9{xi > 1/1) In — 

yi 



[yi + 2/2] - yi 

^^ 1. _^ A (9i-x,)9ix,>y,) + 9ix^)9iy,>x,)) (- + -). (175) 
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Appendix C 



Here, we following Refs. fr7\ ITH] , present the details for the derivation of non-diagonal part of anoma- 
lous dimensions matrix for fermionic operator, based on the analysis of conformal Ward identity for 
Green functions with fermionic operator insertion. In what follows, we will consider somewhat more 
general fermionic operator, then one considered in the main body of the paper: 

ot[ = M^d^ypl'"'^(^)^- (176) 



Even so, this operator will be highest vector in corresponding conformal representation only for a = 1 
and (3 = 0, some of the formulas given here will certainly have applications in other supersymmetric 
models with more rich field content. To start with, let us remained reader some basic facts about 
conformal symmetry group. 

Conformal group is defined as the most general group leaving invariant light-cone in n-dimensional 
Minkowski space. The algebra of this group contains, in addition to generators of Poincare group: 

and Mfj_u, the generator of dilations D and n generators Kfj_ of special conformal transformations. 
The latter generators satisfy the following commutation relations 



[D,K^] = iK^, [K^,P,] = -2t{g^,D + M^,), 

[K^,K,] = 0, [D,M^,] = 0. (177) 

Here we will work only with the subalgebra of conformal algebra, the so called collinear conformal 
algebra, isomorphic to SU{1, 1) = S0{2, 1). The representations of collinear conformal algebra could 
be obtained from those of full conformal algebra via a projection to the light cone. The necessary 
projection itself could be build with the help of two light-cone vectors, n and n* with the properties: 
^2 _ _ g gj^^ n-n* = 1. The generators of collinear conformal group obtained have the following 
expressions: 

K_ = n;K^, D, 

= n^P^, M_+ = nlM^^n^. (178) 

In what follows, we will also need the transformation laws of fermion and scalar fields under special 
conformal transformations: 

K_(j) = i(2x_{d^ + X ■ d) — x'^d_)(f), 

K_il) = i{2x.{d^ + x-d)+2J:^^x'' -x^d.)^p (179) 

where dcf, and are scale dimensions for the scalar and fermion fields correspondingly. is the 
spin operator: 

Sm-^ = ^[7^,7.]^. (180) 

The fermion conformal operators, we are interested in, form an infinite dimensional representation 
of collinear group. The members of each series are labeled by their spin /. The generators D and 
M_+ are diagonal with respect to representations, whereas P+ acts as a raising operator and as 
a lowering operator. The action of on O^^^ is 

K.O-l = / i^. ^ + (KM.)) ^} Ofr (181) 
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Now, it is an easy exercise to find that 

K_Oj7 = a(j,/,a,/5)Oj7_i, (182) 

where 

a{j, I, a, = 2(j - + / + « + /3 + 1). (183) 

To derive the conformal Ward identity itself we start with generating functional for renormahzed 
disconnected Green functions with operator insertions 

Z,MV,J) = ^ y"/^$[Oj7]-exp|2[S]+z Jd''x[vix)i;ix)+ij{x)v{x) + J{x)<P{x)]Y 

i[S] = i I d'^xC, (184) 



where [O^^l^] and [S] denote corresponding renormahzed quantities. Using the invariance of the gen- 
erating functional with respect to special conformal transformations and performing differentiation 
over field sources, one gets the following conformal Ward identity 

([OJT] (^""x)) = -(('^^[OSn) x) - ([OjT] {S''[S]) X), (185) 

where (A) stands for vacuum averaging of the time ordered product TAexp{i[S]) and x = ^i4'i 
denotes product of elementary fields taken at different space-time points. Note, that the left-hand 
side of Eq. p85p is finite, so the right-hand side also should be finite. The first term at right-hand 
side is a variation of the renormahzed fermion operator under special conformal transformations. It 
is given by 



^^[OlT] =zJ2{zMl,e)Z~'}.jOt;U, (186) 

fc=0 "'' 

where we introduced the following notation ^ 

[Otn = i2zA^^9)Ofj, k<l, (187) 



fc=0 



au' 



d^ 



with 



A{1,€) = a(/) + 2(7$ + e)6(0, (188) 
d{l) = {aij,l,a,P)5,k}, (189) 



£ = l-x-^, w(x\a,(3) = (l-xr(l + xf, 
dx 

AfdaB) = 2"+^+^ r(fc + « + l)r(fc + /3 + l) 

i2k + a + P + l)T{k + l)T{k + a + P + l)- ^ ' 



^Thc rcnormalization of conformal operators does not depend on I due to Lorentz invariance and the matrix Zjk 
is triangular. 

^Note, that here we used the fact that anomalous dimensions of scalar and fermion fields are equal, as they belong 
to the same supermultiplet. 
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The evaluation of integral in the definition of matrix b{l) gives 

- ^0>^)|(2^ + « + /^ + l)r(, + « + ;3+l) Tik + a + l)^nk + p + l)\j- 

The second term on the right-hand side of Eq. ()185|1 is the product of two renormalized operators: 
[OjY] and the variation of renormalized action. For the latter we have the following expression 

■^K^g^ = ^iMl[A_] with [A_] = g^t j (fx2x.C{x), ^{e,g) = -eg + (3{g). (192) 

As it is known, the product of two renormalized operators may still contain ultraviolet divergences, 
so an additional subtraction of divergences is necessary 

[0JTA_] = [0J7][A_] - z J2 ZUmOt'J-i]- (193) 



fc=0 



Substituting Eqs. ()18fi|l . ()192j) and ()193|1 into conformal Ward identity and taking limit e 0, one 
gets^ 

{[O^nS-X) =^j2i^{l)+n9,l)hk{[OfJ]x) - ^([0J7A_]). (194) 

fc=0 ^ 

Here the definition for the matrix of special conformal anomalies was introduced 

Y = lim (^ZAZ-^ -a- ^Z*^ . (195) 

It is easy to derive, that in Wess-Zumino model we are considering here Z*^^\l) (the coefficient of 
simple pole in e) is given by 

Z*m^l)=g^(%^l) + 0{g') (196) 

and as a consequence we have Y = bil)i¥^^ + 27$ ^i) with i = {5jk}. 

Now, the crucial point is to note, that the scale and special conformal anomalies are in fact 
related. The relation between matrices of scale and special conformal anomalies could be obtained by 
considering the action of the commutator of special conformal transformation with dilation [6^, 6^] = 
6^ on Green function with fermion operator insertion. Taking into account that 6^ acts on {[0^^l^]x) 

' Oil 



as + iV^7^ + A'"()<,70 (A^^ and A^^ are numbers of i^) and fields in x) one gets the following identity 



«(o+mo+2^KO,7(^?) 

9 



0. (197) 



Since aiV) does not depend on coupling constant, solving this equation recursively we see that 

in (n)-loop order is determined by the (n — l)-loop approximation of 7^^, /5 and 7 = 7^+ 7^^. At 

two first order we have 

,{0) _ r 



7].^ = 5,,7.", (198) 

(0) (0) 

2{j-k)ij + k + 2) 



K.'^r - ... \^I\^.. (^r-/^oM- (199) 



^Note, that at this point the finiteness of the right-hand side of CWI is used 



37 



